Noname manuscript No. 

(will be inserted by the editor) 



Upper and lower bounds for an eigenvalue 
associated with a positive eigenvector 

£>"y Amaury Mouchet 

^3 '■ mouchet@phys.univ-tours.fr 

Laboratoire de Mathematiques et de Physique Theorique (CNRS umr 6083), Uni- 
versite Francois Rabelais Avenue Monge, Pare de Grandmont 37200 Tours, France. 

<n ; 

February 1, 2008 

Abstract When an eigenvector of a semi-bounded operator is positive, we 
show that a remarkably simple argument allows to obtain upper and lower 
bounds for its associated eigenvalue. This theorem is a substantial gener- 
alization of Barta-like inequalities and can be applied to non-necessarily 
purely quadratic Hamiltonians. An application for a magnetic Hamiltonian 
is given and the case of a discrete Schrodinger operator is also discussed. It is 
shown how this approach leads to some explicit bounds on the ground-state 
energy of a system made of an arbitrary number of attractive Coulombian 
i particles. 
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1 Introduction 



In most situations, the principal eigenvalue of a semi-bounded operator 
cannot be obtained explicitly whereas it plays a crucial role in physics: the 
smallest vibration frequency of an elastic system, the fundamental mode 
of an electromagnetic cavity, the ground-state energy of a quantum system 
with a finite number of degrees of freedom, the energy of the vacuum in a 
quantum field theory, the equilibrium state at zero temperature in statistical 
physics, etc. There are actually very few ways — which are usually specific 
to a restricted class of systems 1 — to obtain accurate approximations of an 
eigenvalue with a rigorous control on the errors and a reasonable amount 
of numerical computations. For instance, in a typical Dirichlet-Laplacian 



1 Finding lower bounds for the smallest eigenvalue of a typical Hamiltonian is far 
much difficult than finding upper bounds. For successful attempts, see for instance 
the moment method proposed in |HB85| . the Riccati-Pade meth od propose d in 
FMT89. and the lower bounds obtained for few-body systems in BMB+98 . 
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problem defined for an open connected set Q C K d ; d £ N; Barta's inequali- 
ties Bar37 allow to bound to the lowest eigenvalue eo: The determination of 
a lower (resp. upper) bound requires the finding of the absolute minimum 
(resp. maximum) of a smooth function defined on Q. Compared to the 
general and traditional methods like the Ray lcigh- Schrodinger perturbative 
series and the Rayleigh-Ritz or Temple variational methods, an advantage 
of Barta's approach is not only to naturally provide both an upper and a 
lower bound, but also does not involve the calculation of any integral. There- 
fore, generalizations of Barta's inequalities can lead to interesting spectral 
information. This generalization has been carried out in two directions: 

(i) for Laplacian operators acting on square integrable functions defined 
on a Riemannian manifold (for a recent work on this subject, sec BM04 ); 

(ii) for Schrodinger operators of the form — A + V acting on square in- 
tegrable functions defined on an open set of W 1 _Bar 78 , Bau79 . Thi79 , CR82 , 
Sch85 or more generally second order elliptic operators PW66.BNS94 

EEU. 

In both cases, the proofs of Barta's inequalities involve the use of a 
Kato-like inequality and therefore rely extensively on the somehow specific 
properties of the purely quadratic differential operator. In this paper, I 
propose a significant extension of Barta's inequalities that will rely on the 
properties of one eigenvector only. More precisely, with a remarkably simple 
argument, I will show that we can obtain upper and lower bounds for the 
eigenvalue eo associated with the eigenvector <Pq of an operator, under the 
only hypothesis that <Pq is real and non-negative. This result includes cases 
(i) and (ii) because the Krein-Rutman theorem guarantees the positivity 
of <Pq RS78a §XIII.12] for the smallest eigenvalue eo 2 , but also applies for 
a much wider class of operators including 

(iii) the Schrodinger operators involving a magnetic field, e.g. the Ry- 
drogen atom in a Zeeman configuration; 

(iv) discrete Bamiltonians e.g. the one occurring in the Barper model; 

(v) some integral or pseudo-differentiable operators, e.g. the Klein-Gordon 
or spinless Salpeter Bamiltonians. 

The next section fixes the notations, proves the main general results 
(theorems ^ and |2J) . Section [21 shows how the original argument given in 
section |21 actually embraces and generalizes the Barta-like inequalities that 
have been already obtained in the literature and furnishes guidelines for 
numerically improve the bounds on eo- Sections 0] and El provide two appli- 
cations in the differentiable case (many-body problem) and in the discrete 
case respectively. 



2 The positivity of <Po is also required for the traditional proofs of Barta's in- 
equalities; this explains why, in case (i) and (ii), they concern the lowest eigenvalue 
only. 
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2 Bounding the principal eigenvalue with the local energy 

2.1 General inequalities 

In the following, Q will be a locally compact space endowed with a posi- 
tive Radon measure fi. (ip\<p) will denote the scalar product between two 
elements ip and <p belonging to the Hilbert space of the square integrable 
complex functions L 2 (Q,,u): 



Q 

D(H) will denote the domain of the operator H acting on L 2 (Q, dfi). 
The crucial hypothesis on H is the following: 

Hypothesis 1 The operator H is symmetric and has one real eigenvec- 
tor <Pq G D(H) such that <Pq ^ (almost everywhere with respect to p.) 
on Q. 

If eo stands for the eigenvalue of H associated with <Po, the symmetry 

of H implies that, for all ip G D(H), we have (<Po\(H — eo)(p) = that is 

Vp G D(H), f ¥^(q) (H - e )<p(q) dfi(q) = . 

Jq 

Taking the real part of the integral, we can see that the support of q i— > 
Re I $o (q) (H — eo)p>(q)\ either is empty, either contains two disjoints open 
sets Q± such that Re((Po(H — e )p) and /x(Q±) > 0. The hypothesis of 
positivity of <P implies that on Q±, we have Re[(iJ — eo)<£>] ^ 0. The last 
results motivates the following definition: 

Definition 1 (Local energy) For any ip in 0(H), the local energy is the 
function E v : Q — > R defined by 



Rc(g^(g)) 

Re(<p(q)) 



EM = 3 ^ • (i) 



Therefore from what precedes, we have obtained the main theorem: 

Theorem 1 For any symmetric operator H on L 2 (Q,/j.) having an eigen- 
value eo whose corresponding eigenfunction is nonnegative almost every- 
where on Q, we have 

yip G D(H) such that Re(cp) > 0, mf(£L) < e < sup(E v ) . (2) 

2 S 



Actually, for a non-symmetric operator K, we can keep working with its 
adjoint K* and easily generalize the above argument: 
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Theorem 2 Let K being an operator on L 2 (Q,/i) having an eigenvalue kg 
whose corresponding eigenfunction is real and nonnegative almost every- 
where on Q, we have V</? G D(K*) such that ip > 0, 



inf 

Q 



Re(K*(p) 



^ Re(fc 



inf 

a 



Im(fco) 



^ sup 
Q 



^ sup 
a 



Re(iTV) 
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(3a) 



(3b) 



This generalization may be of physical relevance. There are some models 
(e.g. the so-called "kicked" systems, or quantized maps) where the dynam- 
ics is described "stroboscobically" i.e. implemented by a unitary operator 
(the Floquet evolution operator) that cannot be constructed from a smooth 
Hamiltonian. However, we will not consider this possibility here, and up to 
the end of this paper, H will denote a symmetric operator. 



2.2 Optimization strategy 

Since generally, the eigenfunction <Pq is not known exactly, it will be ap- 
proximated with the help of test functions that belong to a trial space 
T(H) C D(iJ), very much like the variational method. Since we want a test 
function to mimic <Pq at best, we will restrict T(H) to functions that respect 
the a priori known properties of its positivity, its boundary conditions 
and its symmetries if there are any. For each test function the error on cq 
is controlled by inequalities (J2J- Therefore, the strategy for obtaining rea- 
sonable approximations is clear: First, we must choose or construct ip to 
eliminate all the singularities of the local energy in order to work with a 
bounded function and second, perturb the test function in the neighborhood 
of the absolute minimum (resp. maximum) of the local energy in order to in- 
crease (resp. decrease) its value. For practical and numerical computations, 
this perturbation will be implemented by constructing a diffeomorphism 
A — > T(H), A i— » ip\ from a finite dimensional differentiable real manifold A 
of control parameters A and the optimized bounds for eo will be 

supinf(£' ¥ , A ) < e < inf sup(£^J . (4) 
A Q A Q 



3 Inequalities in the differentiable case: old and new 

3.1 General considerations 

When H is a local differential operator i.e. involves a finite number of 
derivatives in an appropriate representation (for instance in position or in 
momentum representation), one can therefore construct an algorithm that 
does not require any integration but differential calculus only. In an ana- 
lytic and in a numerical perspective, this may be a significant advantage 
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on the perturbative or variational methods even though it is immediate to 
see 3 that the upper bound given by J5J) is always larger than {ip\H<p) / (ip\ip). 
The global analysis appears only through the determination of the singular- 
ities and the absolute extrema of the local energy that may have bifurcated 
when the control parameter A varies smoothly. For a Schrodinger operator, 
the possible singularities of the potential on Q like a Coulombian diver- 
gence or an unbounded behaviour at infinite distances may furnish a strong 
guideline for constructing relevant test functions (see § 13.41 below) . I have 
given m MouOSb , some heuristic and numerical arguments to show how 
this strategy can be fruitful. In the present article, the main focus will con- 
cern rigorous results and I will explain how some of them can be obtained 
with great simplicity even for systems as complex as those involved in the 
many-body problem. 



3.2 Case (i) Barta's inequalities 

They immediately appear as a particular case of theorem ^ 

Theorem 3 (Barta 1937) Let Q be a connected bounded Riemannian ma- 
nifold endowed with the metric g and H , the opposite of the Laplacian A g 
acting on the functions in L 2 (Q,/i) that satisfy Dirichlet boundary condi- 
tions on the boundary dQ. The lowest eigenvalue eo of H is such that, for 
all positive ip 6 C 2 (Q), 



Proof: It follows directly from theorem ^ with the local energy given 
by E v = —Agip/ip: The spectrum of H is discrete and the Krein-Rutman 
theorem assures that hypothesis Q is fulfilled for eo being the lowest (and 
simple) eigenvalue. □ 

Remark 1 The Dirichlet boundary conditions are not essential and can be 
replaced by any other type of boundary conditions provided that hypoth- 
esis Q] remains fulfilled. However, as explained in § 12.21 for obtaining inter- 
esting bounds on e extending T{H) to test functions that do not fulfill the 
boundary conditions (as proposed in |Duf47| ) seems not appropriate. 



3 For each positive tp, it follows from 



(<p\H V ) = Be[(<p\H V >)] = / <p(q)Be[H<p(q)]dii(q) 



<P (<l)E^(q)da(q) < (<p\tp) sup (E v (q)) . 
Q Q 
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3.3 Case (ii) Duffin-Barnsley-Thirring inequalities 

Extensions of Barta's inequalities for Schrodinger operators have been ob- 
tained partially by Duffin |Duf47| and Barnsley BarTgj (for the lower bound 
only) and completely (lower and upper bound) by Thirring |Thi79| using 
Kato's inequalities (see also Sch85 ). 

Theorem 4 (Duffin 1947, Barnsley 1978, Thirring 1979) 

Let H = —A + V be a Schrodinger operator acting on L 2 (lR d ) having an 
eigenvalue below the essential spectrum. Then the lowest eigenvalue eo of H 
is such that for any strictly positive <p G D(H), 

inf [V - — | < e s$ sup [V - — ) . (6) 

» d V f J R d V <P J 

The proof is similar to the one presented above with the local energy being 
now Eu> = V — Aip/ip. This argument has the advantage on the existing ones 
that it does not involve the specific properties of the Laplacian and can be 
immediately transposed to the larger class of the differential operators (not 
necessarily of second-order) that fulfill hypothesis ^ 



3.4 Case (Hi) Magnetic Schrodinger operators 

In the presence of a magnetic field, Schrodinger operators take the form 
H = (id q + A(q)) 2 + V(q) with A : Q — > R d being a smooth magnetic 
potential vector and V : Q — * K a smooth scalar potential. The Krein- 
Rutman theorem may not apply whereas there still exists a non-negative real 
eigenfunction (e may be not simple nor the lowest eigenvalue) HHH099 . 

In the particular case of the Hydrogen atom in a constant and uniform 
magnetic field, hypothesis 03 is fulfilled for all values of the magnetic field 
AHS77 AHS78 and indeed concerns the lowest eigenvalue. Therefore theo- 
remHapplies and furnishes relevant analytical bounds that can be improved 
numerically as shown in [MouORbj . 

Proposition 1 The smallest eigenvalue eo of the (Sd-)Zeeman Hamiltonian 
!^I(-iV + irxB) 2 -i (7) 

is such that 

VB^O, e ^ -1/2 + B/2. (8) 

Proof: In cylindrical coordinates (p, 9, z) where the magnetic field 
is B = Bu z , the test function of the form tp = cxp(—^/p 2 + z 2 — Bp 2 /4) 
is constructed, according to the strategy explained in § 12.21 in order to 
respect the rotational invariance of the ground-state AI1S77 and eliminate 
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both singularities at r — > and at p — > oo. Indeed such a choice leads 
straightforwardly to the bounded local energy: 

E v = -1/2 + B/2 - - (9) 



2v^ 2 + 



The upper bound follows. □ 



4 Application to the many-body problem 

4-.1 Expression of the local energy in terms of two-body functions 

We will consider in this section a iV-body non-relativistic bosonic system in 
d dimensions; (d, N) £ (N\{0, l}) 2 ; whose Hamiltonian is given by: 

N-l 

acting on L 2 (]R Ar<i ), endowed with the canonical Lebesgue measure, and 
where Vi € {0, . . . , N — 1}, r,* S R d . Z\j is the Laplacian in the variables 
and mi € R + \{0} the mass of the i th particle. The spinless bosons interact 
only by the two-body radial potentials Vij = Vji : R — > K i.e. V is given by 

N-l 

i<j 

where = rjj = ||rj — r$||. Once the center of mass is removed, the 
Hamiltonian H leads to a reduced Hamiltonian H acting on L 2 (R( Ar_1 ' ci ) 
(see for instance § XI.5 of RS78b ) and we will suppose in the following that 
H has at least one eigenvector 4 . Therefore, hypothesis Q is fulfilled for eo 
being the lowest (and simple) eigenvalue. A natural choice for test functions 
is to consider factorized ones of the form: 

N-l 

<p(r ,...,T N -i) = Yl 4>ij(nj) (12) 

i,3=0 
i<j 

4 Physically, this can be achieved with a confining external potential (a "trap" 
is currently used in experiments involving cold atoms). Formally, this can be 
obtained in the limit of one mass, say mo, being much larger than the others. The 
external potential appears to be the «o,i's, created by such an infinitely massive 
motionless device. It will trap the remaining N — 1 particles in some bounded 
states if the uo,i's increase sufficiently rapidly with the ro.i's. 
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where <pij £ L 2 (R + ) and cj>ji — (f>ij > 0. One can check easily that the 
total momentum of such a test function vanishes. The corresponding local 
energy is given by 



N! 



N-l 

E 

i,j=0 



4>ij ( r *i ) + - — - 4>'ij {n, ) ) + v i] (n. 



2mij(f>ij(rij) 

i<3 (13) 
~ E — S 'ij( r ij) S ik( r ik)cos(j,i,k) 

where qjv € Qn = R^ -1 ^ stands for the (N — l)d relative coordinates 
(rj. — ro, . . . , r n— i ~ro), ffiy for the reduced masses mirrij / (rrii + rrij) , S'^ is 
the derivative of Sij = ln(0y). The last sum involves all the N(N — l)(N — 

2)/2 angles (j, i, k) between Tj — and r& — that can be formed with all 
the triangles made of three distinct particles. 

Whenever each Vij allows for a two-body bounded state, we can choose 
4>ij to be the eigenvector of — (2m !3 -)" 1 Zl + having the smallest eigen- 

(2) 

value ■ . Moreover, if Vij(r) is bounded when r — > oo, from an elementary 
semiclassical analysis (see for instance [MES11 ) it follows that S[j is also 

bounded since asymptotically we have (r) ^ r ->oo — \J 2rnij [v^ (r) — ej^] . 
It follows that the local energy is also bounded and finite lower and upper 
bounds on eo can be found. For instance, directly from expression Q13[l. we 
have 



Proposition 2 //, for all £ {0, .. ,,N — l} 2 , i ^ j, Vij(r) is bounded 

Vij has a smallest eigenvalue e^ 1 



when r — > oo and -(2mjj) 1 A + Vij has a smallest eigenvalue obtained 



for 4>ij — expSij, then the smallest eigenvalue eo of the N-body Hamiltonian 
in the center-of-mass frame is bounded by 

N-l 2 N-l 2 

E i 2) -^^-l)(^-2) < eo < £ 4 2) + ^7V(iV-l)(7V-2). (14) 

iJ—O 

where m = min; mj and s = max^j sup R + |Sj,-|. 

For potentials that are relevant in physics (see for instance the effec- 
tive power-law potentials of the form v g(r) — sign(/3)r^; (3 £ R; between 
massive quarks as studied in |BMB + 98| '). the analytic form of the two-body 
eigenvector is not known in general and some numerical computations are 
required to obtain the absolute maximum and minimum of the local en- 
ergy (|T3)l . 



4-2 The local energy for a general Coulombian problem 

When N is large, the estimation l|14f> is quite rough, in particular it does not 
take into account the constraints between the several angles. More precise 
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results are obtained for the Coulombian problem where Vij(r) = Cij/r with 
eij G M. and d > 1. In that case, provided a bounded state exists and 
that we keep the test function 1)12(1 in L 2 (R^ Ar_1 ^ d ), we choose a constant 
derivative SL = 2eijrriij / (d — 1) in order to get rid of the Coulombian 
singularities of V . We obtain a bounded local energy given by 

Mm) = E q - (rf31)2 - (jTija £ ^ «*(,,.,*). (15) 



4-3 Identical purely attractive Coulombian particles 

The case where all the N particles are identical and attract each other, i.e. 
when V(i,j) S {0, . . . , N — l} 2 , i ^ j, mi — 1 and = —1, has been 
extensively studied in the literature, in particular the asymptotic behavior 
of eo with large N may have some dramatic consequences on the thermody- 
namical limit FR6 6|ILD67IILL69| . The local energy method allow to obtain 
in a much simpler way energy bounds that are comparable to those already 
obtained by other methods. Actually, (|15|) simplifies to 



with F 2 = and for N > 3 



^N{N-l)+F N (q N ) 



(16) 



Fn(qn) = cos U,i,k)- (17) 

(j,i,k) 

The angular function Fn depends only on the geometrical configuration of 
the N vertices (ro, . . . , rjy-i) i.e. it is invariant under the group of Euclidean 
isometries and the scale invariance of the Coulombian interaction makes it 
invariant under dilations as well. 

Lemma 1 inf q 3 F3 = 1 is obtained when the three points are aligned. supg 3 F3 
3/2 is obtained when the three points make an equilateral triangle. 

Proof: The extrema of F3 correspond to the extrema of the function defined 
by (#oj 6*i , 6* 2 ) 1— > cos 8q + cos Q\ + cos 9 2 under the constraint 8 + 61 + 8 2 = 7r 
for (0 , 81,82) being the three angles (1,0, 2), (0, 1, 2), (0, 2, 1) respectively. 
The Lagrange multiplier method leads to the determination of the extrema 
of the function [0; 7r] 3 — > R defined by (80, 81,82) ^ cos 9q + cos 8\ + cos 82 + 
1(80 + 81+82 — it); I € K. We immediately obtain that the extremal points are 
located at (80,61,82) = (tt/3, ir/3, n/3) and (80,81,82) = (0,0, n) together 
with the solutions that are obtained by circular permutations. It is easy to 
check that the first solution provides an absolute maximum for F3 and the 
second ones an absolute minimum. □ 

An immediate consequence of the preceding lemma is 
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Proposition 3 The lowest energy eo of N ^ 2 identical attractive Coulom- 
bian spinless particles in d > 1 dimensions is such that 

e ^--(jh^l N{N - 1){N+1) (is) 

w/iera f/ie individual masses equal to unity and the attractive potential is — 1/r. 

Proof: The sum on the angles that defines Fn in (|17fl can be written as a sum 
of N(N — 1)(N — 2)/6 i^-terms calculated for all the triangles that belong 
to the iV-uplet made of the N vertices. Since, from Lemma Q] the absolute 
minimum of F3 is obtained for a flat configuration, when all the N points 
are aligned all the i^-terms reach their absolute minimum simultaneously 
and the absolute minimum of Fpf is obtained. We have 

Lemma 2 infq^ Fjf = — N(N— l)(N— 2) is obtained when all the N points 
are aligned. 

The upper bound (|18|) follows from i|16|) . □ 

More generally, for a given 7V-uplet (i.e. a set of exactly N points), 
clustering the sum JT7j > in M-uplets (N > M ^ 3) allows to find bounds 
on Fn from bounds on Fm- Indeed, we can write 

, . (M — 3)! (N — M)\ , , 

F N (q N ) = ]T ^~3M ~ F m(<1m) (19) 

qM ^ '' 

where the sum is taken on all the M-uplets, labeled by the coordinates qM, 
that belong to the given iV-uplet. This sum involves exactly N\/M\/(N — 
M)\ terms and we have, therefore, 

Lemma 3 V(iV, M) € N 2 such that N > M ^ 3, 

N(N-l)(N-2) , s 

For a given N, supq n F^ is not known exactly but the ordered sequence 
sup Cjv F N < SU PQ W -! f n-i < 



N(N - 1)(JV - 2) (iV - l)(iV - 2) (AT - 3) 

< su PQm F m su Pq 3 f 3 = 1 ( 

^ M(M- 1)(M- 2) ^ ^ 3.2.1 4 1 j 

shows that in order to improve the lower bounds on (|16|l , we must try to find 
supq flf Fm with M being the largest as possible. However, when considering 
identity (|19fl for N = 4 and M = 3 together with Lemma ^ we have 



Lemma 4 supq 4 -F4 = 6 is obtained when the four points make a regular 
tetrahedron. 
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Then supg 4 i<4/(4.3.2) = supg -^3/(3. 2.1) = 1/4 and no better estimate is 
obtained when considering M = 4 rather than M = 3. Numerical investi- 
gations lead to the following conjectures: 

Conjecture 1 (C$) When d = 3, supF 5 — ^ + = = — -— ^ — j- 



2 + § 



14.591594 with 



6/i = l + W-l+ V / 7 + 4v / 3 



V7 + 4V3 



7 + 4^3- 



^ ^7 + 4^3 ^/-l + (7 + 4 V3)V3 + (7 + 4 V3)- 1 / 3 

(22) 

is obtained when the five points make two mirror-symmetric tetrahedrons 
sharing one common equilateral basis, their other faces being 6 isosceles 
identical triangles. 

Remark 2 The only free parameter of the specific configuration can be cho- 
sen to be the height h of one tetrahedron (the length of the edges of the 
common equilateral basis being fixed to one). The maximum of h 1— ► is 
reached for ho being the greatest solution of 9 ft 4 — 6h 3 + 3h 2 — 2h + 1/3 
that is precisely given by (|22|) . 

Remark 3 The pyramidal configurations with a squared basis leads to a local 
maximum that gives F5 = 15/2 + 5\/2 ~ 14.57. 

Conjecture 2 (Cq) When d — 3, supg g F& = 12(1 + v2) is obtained when 
the 6 points make a regular octahedron. 

Conjecture 3 (C s ) When d = 3, su Po F 8 = ie[| + -4 + -4 + 4 ( 1+ ^) + 

^ s V 2 v5 V5V5+4V2 

3+2y/2 ^ 79 5QX i s obtained when the 8 points make two identical 

squares (whose edges have length one) lying in two parallel planes separated 
by a distance h = yl + 2 v / 2/2. The axis joining the centers of the two 
squares is perpendicular to the squares and the two squares are twisted one 
from the other by a relative angle of 7r/4. 

Remark 4 The cube corresponds to F s = 8(3^2 + a/3 + 3/2 + ~ 79.393. 

Conjecture 4 (Coo) When d = 3 and AT — > 00, the configuration that max- 
imizes Fn corresponds to N points uniformly distributed on a sphere and 
su PQn F n ~IN 3 + o(N 3 ). 

Remark 5 The ambiguity of distributing N points uniformly on a sphere 
[SK97I and references therein] vanishes for large N as far as a uniform 
density is obtained. Assuming such a uniform density, the continuous limit 
of Fn/N 3 is a triple integral on the sphere than can be computed exactly 
to 2/9. 
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The upper bound IjlSjl for d = 3 is slightly above the one obtained 
in |LL69I eq. (17 p. 807)], na mely -(5 /8) 2 iV(jV - l) 2 /8, and the numerical 
estimate -.05427V(iV-l) 2 in BMR90, eq. (16) p. 63]. For the lower bounds, 
from Lemmas n or we have obtained 

Proposition 4 The lowest energy eo of N ^ 2 identical attractive Coulom- 
bian spinless particles in d > 1 dimensions is such that 

' ]n 2 (N-1) ^e (23) 



(d-l) 2 4 

when the individual masses equal to unity and the attractive potential is — 1/r. 

The same result has been obtained for d = 3 in BMR90, eq. (11) p. 62] 
and is slightly better than -N(N - l) 2 /8 given in |LL691 eq. (13 p. 807)]. 
For N ^ M ^ 4, this lower bound can be improved to 

1 N(N-l)(l+a M (N-2)) ^e Q (24) 



[d-lf v ,y 2 
with 

sup Qju F M 1 . . 

aM ~ M(M-l)(M-2) ^ 4- (25) 

If conjecture (C5) [resp. (C*6), (Cs) and (Coo) ] is correct we get, when d = 3, 
as ~ .2432 [resp. a 6 ~ .2414, a 8 ~ .2366 and ~ 2/9 ] and the lower 
bounds are, therefore, improved. 

Some numerical investigations, in particular a systematic comparison 
with the lower bounds obtained with variational methods in |BMB + 98] for 
N = 3 and N — 4 Coulombian particles will be given elsewhere jMou05a . 



5 Application to discrete Hamiltonians 

Where q £ Z d ; d € N; the discretized analogue of a local differential operator 
corresponds to a Hamiltonian that couples at most a finite number of basis 
vectors (e.g. the nearest neighbors on the lattice 1 d ). For instance, when d = 
1, it can be seen as a Hermitian band matrix (finite or infinite) of finite half- 
width in an appropriate basis. Possibly with renumbering the q's, on £ 2 (Z d ) 
H has the form: 

i Ht P)q = Hq-q+v^q+v > ( 26 ) 

where N b G N, V(g,q') e Z d x Z d , = 5,,^ € C and stands for 

max(|^i|, . . . , l^rfl). <p will be taken as a discrete set of real strictly positive 
numbers, and the local energy E v (q) is computed, for a given q, with el- 
ementary algebraic operations whose number is finite and all the smaller 
than N b is small: its value at a given q depends on (2N b ) d + 1 components 
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of ip at most. Under hypothesis^, if, say, for a given test vector <p, the 
absolute maximum of E v occurs only at a unique finite q m , one can imme- 
diately improve the upper bound by a finite amount, for instance just by 
varying (p qm only, until E v (q m ) is not an absolute maximum anymore. Only 
(2N b ) d + 1 values of the local energy will be affected by the variation of just 
one component of ip. One can see easily that this approach leads to a wide 
variety of algorithms where a sequence of optimization steps is constructed; 
each step involves a number of optimization parameters and functions that 
is usually much smaller (of order (2Nb) d or less) than the dimension of the 
original matrix. 

Discrete Schrddinger operators are important particular cases of Hamil- 
tonians (|26|) with Nf, = 1. They are relevant models for the description of 
quantum (quasi-)particles evolving in periodic crystals. For d = 1, they can 
be written as 

(H<p) q = -tp q +i - tp q -i + V(q) <p q (27) 

where the potential V is a real bounded function on Z. By possibly subtract- 
ing a constant positive real number to V, the bounded operator — H can be 
made positive and ergodic: First, for any positive and non-identically vanish- 
ing ip and ip' in ^ 2 (Z), —Hip remains positive. Second, H)\ q ~ q][ ip) = 
f'q'Vq + (positive terms) ^ for any given pair of strictly positive com- 
ponents ip' q ' and ip q with q' ^ q (In the marginal case where ip and ip' 
both vanish everywhere but on the same point, we have (ip'\(—H)ip)) ^ 0). 
Therefore theorem XIII. 43 of RS78a applies: if H has indeed one eigen- 
value, hypothesis n is fulfilled for eo being the smallest eigenvalue of H. In 
that case, inequalities (J2J) takes the form: 

Proposition 5 When the discrete Schrddinger operator l|27(l admits at least 
one eigenvalue and when V is bounded, then the smallest eigenvalue eo is 
such that, Vip £ £ 2 (Z) such that (p > 0, 



inf ( - 



+ V(q)) < e < sup f- jg*H±jgg^ + V(q)) . (28) 




When V is actually a iV-periodic real function, the spectral problem (see 
|RS78a| for instance) leads to the search of complex series (u q ) q ^z such that 



ir , - . -Uq- 1 + V{q + m)Uq = e{rf)u q ; 
Vry=(r ? i,r ?2 )G[0;l[ { (29) 

The spectrum of H is the bounded set a(H) = {e(n)\i] G [0; 1[ 2 } C K. 
It is given by the reunion for all ?7's of the N eigenvalues of finite N x N 
Hermitian matrices obtained after transforming (|29|l with the one-to- 
one mapping u q i— > u q exp(—i2Trqrii/N). As far as positive solutions of l|29|l 

5 We have seen at the end of section \2. II that the symmetry hypothesis can be 
relaxed. 
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are concerned, we will take 771 
value eo(r)2) of 

fv(m) -1 

-1 V"(l + 772) -1 • 
ff (0 <' ?2) =0-1 V(2 + r] 2 ) -1 . 



\ -1 ••• 0-1 V(N-l + m )J 



and will look for the smallest eigen- 
■•• -1 \ 

(30) 



-1 




Remark 6 The (rational) Harper model Har-55b,Har55a (also called the 
almost Mathieu equation BS82 ) corresponds to V(q) = — Vq cos(2irqM/N) 
where Vq > 0, (M, N) being strictly positive coprimes integers. For a given 
N and M, cr(H) appears to be made of N bands. The union of these bands 
for each rational number M/N between and 1 produces the so-called Hof- 
stadter butterfly |Hof76| . 

We are therefore able to produce two non-trivial bounds on the lowest eigen- 
value eo(r]2) without any diagonalization: 

Proposition 6 When V is N -periodic, V772 € [0; 1[, the smallest eigen- 
value eo{r] 2 ) of (EOJ) is such thatVtp € {M+\{0}) N , 



mm 

q e{o,...,N-i} 



Vq+1 + Vq-l 



Vq 



+ V(q + T] 2 ) «S e {r] 2 ) 



en (779) ^ max 

ge{0,...,JV-l} 



ip q+1 + (Pq-1 
Vq 



+ V{q + 772) 



(31a) 
(31b) 



(the indices labeling the components of tp are taken modulo N). 

Therefore, we can bound the bottom of the Hofstadter butterfly with 
the help of any test function: 

Corollary 1 For the rational Harper Hamiltonian 



M 

Htp q = -<p q+ i - ip q -! - V Q cos ( 27rg— 



Vq 



(32) 



we have Mip G (M+\{0}) iv 

Vq+1 + Vq-l 



mm 

ge{o,...,JV-i} 



Vq 



Vq cos 2nq 



M 

'iV 



inf [a{H)} (33a) 



and 



inf \a(H)] ^ max 

g e{o,...,iV-i} 



Vq+1 + Vq-l 
Vq 



V COS [ 27T<7— I 



(33b) 



Proof: In the Harper model, for each rational number M/N, the lowest 
eigenvalue is obtained for 7/ = (0,0). It is a direct application of |TT5 78a 
Theorem XIII. 89 (e)] and thus inequalities (|33|l follow directly from 1)31 [I. 
□ 

Choosing for tp, at first guess, semiclassical approximations (i.e. cor- 
responding to large N) constructed from Mathieu functions are therefore 
expected to provide numerical reasonable bounds. 
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6 Conclusion 

It has been shown on various examples how theorem ^ can be used to ob- 
tain rigorous estimates on the principal value of any symmetric operator. Its 
simplicity, its low cost in computations and its wide domain of applications 
make the method presented in this article a powerful tool for controlling 
bounds. In many situations, it provides non-trivial complementary informa- 
tion to those obtained by traditional or more system-dependent methods. 
Unfortunately, I was not able to extend the method to fermionic systems 
(see for instance |Sig95| and references therein) where the spatial wave- 
function of the ground-state has generically non-trivial nodes |Cep91| that 
cannot be known a priori even with some considerations on symmetries. 

I have also presented some clues for further developments of optimization 
algorithms. However it remains an open question whether such algorithms 
really bear the potential of an efficient treatment and will overcome the 
possible difficulties one may face in realistic problems. 
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